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Q ^ The effects of entanglement in solutions and melts of unknotted ring polymers have been addressed 

by several theoretical and numerical studies. The system properties have been typically profiled as 
a function of ring contour length at fixed solution density. Here, we use a different approach to 
00 ' investigate numerically the equilibrium and kinetic properties of solutions of model ring polymers. 

^Nj , Specifically, the ring contour length is maintained fixed, while the interplay of inter- and intra-chain 

entanglement is modulated by varying both solution density (from infinite dilution up to « 40% 
volume occupancy) and ring topology (by considering unknotted and trefoil-knotted chains). The 
^-H , equilibrium metric properties of rings with either topology are found to be only weakly affected by 

O . the increase of solution density. Even at the highest density, the average ring size, shape anisotropy 

and length of the knotted region differ at most by 40% from those of isolated rings. Conversely, 
kinetics are strongly affected by the degree of inter-chain entanglement: for both unknots and trefoils 
the characteristic times of ring size relaxation, reorientation and diffusion change by one order 
of magnitude across the considered range of concentrations. Yet, significant topology-dependent 
differences in kinetics are observed only for very dilute solutions (much below the ring overlap 
threshold). For knotted rings, the slowest kinetic process is found to correspond to the diffusion of 
^ ' the knotted region along the ring backbone. 
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> ■ I. INTRODUCTION 

m 

' Characterizing the equilibrium and kinetic properties of semi-dilute solutions of cyclic polymers (rings) is one of 
\^ I the major remaining challenges in theoretical and experimental polymer physics [J-Q]- 

One aspect of these systems that is very actively investigated regards how intra-chain and intcr-chain entanglement 
reverberate on physical properties of rings solutions. Entanglement effects have been extensively investigated, and are 
hence well understood, for linear, i.e. open, polymer chains. However, unlike linear chains, physical properties of dense 
solutions of rings strongly depend on their preparation history during which the ring's geometrical entanglement is 
locked permanently in topological constraints such as links and knots. Because of these constraints, the same concepts 
that proved so useful in the characterization of dense solution of linear chains are usually not easily, nor profitably, 
transferred to the case of circular polymers. 
^ ' A chief example, that paradigmatically highlights the qualitative differences of systems of closed versus open chains, 
is offered by the so-called Edwards-DeGenncs reptation model [8|,iS] for the kinetics of dense linear polymers. According 
to this model, each chain is confined inside a tube-like region resulting from the excluded- volume interactions with the 
neighboring chains. The motion of each chain therefore consists of a one-dimensional diffusion along the tube centerline 
resulting from the "inchworm movements" of small sub-chain loops. This type of local motion is found also in polymer 
rings moving through a fixed array of obstacles However, the asymptotic standard diffusive behavior of polymer 
chains in melt depends on the ongoing process of "tube renewal" which results from the fact that the two polymer 
ends are practically free and hence capable of exploring and realizing new configurations This mechanism, 

which has received striking experimental confirmations as being the dominant kinetic process in dense polymer melts 
[l^ . explicitly builds on the linear character of the chains in solution. Hence, it is not applicable, at least in its 
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conventional formulation, to ring, polymers. The latter must consequently relax and move in dense solution according 
to different kinetic mechanisms [2|, which are only recently being characterized by means of computationally-intensive 
molecular dynamics simulations [l, 0, IT5l - [l9l | . 

Ring topology affects strongly also equilibrium properties of the chains, and in particular metric ones (20l - [23 |. In 
fact, the closure condition not only introduces correlations in the backbone orientation that are of longer range than 
for linear chains HI] , but it preserves the knotted topology of the ring. If chain self-crossing is disallowed (no bond 
is cut and rejoined) then the ring knotted or unknotted state, trapped at the time of circularization, is maintained at 
all times. 

The ring self-entanglement associated to a given topological state is so important that it affects the scaling of metric 
properties of isolated rings, as first argued by Des Cloizeaux [23|, and later by Grosberg and Deguchi [2^,[29| and other 
studies (see cited references in [l^l). In particular. Gates and Deutsch and Miiller et al. 0,131 have suggested that 
strong inter-chain topological interactions lead to pronounced chain compaction of unlinked and unknotted rings in 
semi-dilute solutions, and that the typical chain configuration should resemble a branched polymer (akin to a lattice 
animal) [l,!]. 

These seminal studies spurred several theoretical and computational efforts where equilibrium and kinetic properties 
of ring solutions were studied at fixed density - almost invariably corresponding to melt conditions - and for increasing 
ring length. These approaches, have provided considerable insight into the expected asymptotic scaling of the ring's 
metric properties 0, 0, [H, El]- In particular, recent evidence d, 0] suggests that the branched polymer regime is 
likely at a crossover between the Gaussian regime for short rings and the "crumpled globule" regime for long ones 



31|. In the latter regime, each ring portion is highly compact and minimally intermingling with the rest of the chain 
31 1, similarly to dense systems of long biopolymers in vivo, such as eukaryotic chromosomes j32| . 

Gompared to these studies, here we use stochastic molecular dynamics simulations to investigate complementary 
aspects of the effects that intra- and inter-chain topological constraints have on kinetic and equilibrium properties 
of solutions of model ring polymers. Specifically, we consider systems in which the contour length of the rings is 
fixed while the solution density is varied so to cover an appreciable range of inter-chain entanglement. The study is 
carried out for monodisperse solutions of unknotted rings as well as trefoil-knotted ones. To have a better insight 
into topological effects on the statics and dynamics of knotted rings we look also at the knotted portion of the ring. 
This is in general a formidable task and we succeeded in tackling it by relying on a recently-introduced effective and 
transparent procedure for locating knots [H, HJ] ■ To the best of our knowledge this is the first time that a systematic 
off-lattice study of kinetic and equilibrium properties of rings with non-trivial topology and of their knotted portions 
has been carried out. 

By first exploring the metric and shape properties of equilibrated rings in solution wc find that their behavior 
depends appreciably on ring topology and that not only it differs from the one seen for linear polymers in solutions 
but even from other physical realizations of dense ring polymers. Specifically, unlike what has been observed in 
confined f33l - l35| and collapsed knotted rings [H, H^] , the average size of the knotted portion of a ring in solution is 
only weakly affected by the increase (for increasing monomer concentration) of the system's geometrical entanglement. 
In particular, by increasing the solution density, no crossover to delocalizcd knots and no multiscale behavior of the 
entanglement is observed, unlike the case of spherical confinement [s^ . 

The presence of topological constraints makes the dynamics of rings in solution an even richer phenomenon whose 
understanding requires a full knowledge of the relationship between the spatial motion of a ring in a melt and the 
motion of its knotted portion along the ring backbone. By using the knot location tool previously mentioned, and 
considering dynamic observables commonly used in linear polymer contexts we show that several time scales are at 
play: the autocorrelation time of ring configurations, the time required by the rings to diffuse over regions comparable 
to their average size and - for knotted rings - the time required by the knot to diffuse over the ring contour. 

Various relaxation properties of the entire ring are found to change by one order of magnitude across the considered 
density range. The slowest kinetic process is associated with the diffusion along the backbone of the knotted region. 
This property, that to the best of our knowledge has not been pointed out before, seamlessly integrates with the other 
kinetic aspects thus offering a consistent picture for polymer relaxation in solutions of topologically constrained rings. 



II. MODEL AND METHODS 



A. The model 



To model the rings in solution, we use the bead-spring polymer model introduced by Kremer and Grest [38| . The 
model accounts for the connectivity, bending rigidity, excluded volume and topology conservation of polymer chains. 
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Specifically, the intra-chain energy consists of the following terms: 

n 

'Hintra=^[UFENE{i,i + i)+ (1) 
1=1 

Ubr{i,i + lJ + 2)+ (2) 

n 

E ULA^,J)] (3) 

j=i+i 

where n is the total number of beads per ring, and i and j run over the indices of the beads. The latter are assumed 
to be numbered consecutively along the ring from one chosen reference monomer. The modulo-71 indexing is implicitly 
assumed because of the ring periodicity. 

From now on we shall take the nominal bead diameter, a, as the unit length and adopt the following notation: the 
position of the center of the ith beads is indicated by r*; while the pairwise vector distance of beads i and j is denoted 
as dij ~ fj — fi and its norm simply as di,j. 

With this notation the chain connectivity term, UpENEihi + 1) is expressed as: 



UFENE{^,^ + l)= { 2^°^'' 



1 - 

0, rfi.i+i > Ro 



Ro 



where Rq = 1.5cr, k = SO.Oe/cr^ and the thermal energy ksT equals l.Oe [Sl]. The bending energy has instead the 
standard Kratky-Porod form (discretizcd worm-like chain): 

(7f,r(z,i + 1,1 + 2) = 1 ^ ~3 3 (5) 

C \ <J'i,i+l «i+l,i+2 I 

where = 4.5ct is the nominal persistence length of the chain [39l |. Polymer chains are significantly bent by thermal 
fluctuations at contour lengths larger than where = 2^p = 9. Oct is the Kuhn length of the chain d, H^. 

The excluded volume interaction between distinct beads (including consecutive ones) corresponds to a purely 
repulsive Lennard- Jones potential: 

^^■^^*'^^-\0,d.,, >ct21/6 ■ (6) 

This repulsive interaction controls the inter-chain excluded volume too: 

^^^nter = E ^Lj{i,j) (7) 

1=1 ,7=7+1 

where is the number of rings in solution and the index i [j] runs over the beads in chain / [J]. 



B. Simulation details 



We consider solutions of = 64 rings, each consisting of n = 216 beads (i.e., the ring contour length, Lc = 216ct 
corresponds to 2Mk) at six different monomer densities, p: pa^ = 0.010,0.025,0.050 < p*(7^ ■, pcr^ = 0.100 and pa^ = 
0.200,0.400 > p*cr^. Notice, that the chosen densities cover the cross-over to the overlap monomer concentration, 
p*a^ « 0.1 (see Eq. [TUl below) . For reference, we also consider isolated rings, and solutions of = 21 rings in very 
dilute conditions {pa^ = 0.003). 

The equilibrium and kinetic properties of these systems are studied using fixed-volume and constant-temperature 
Molecular Dynamics (MD) simulations. The system dynamics is integrated with the LAMMPS engine [40] with 
Langevin thermostat (target temperature = 1.0 LJ-units). The elementary integration time step is chosen equal to 
At = 0.012rj\775, where tmd = a{rnj_£)^l'^ is the Lennard- Jones time, m is the bead mass and the friction coefficient, 
7, corresponds to 7/m = 0.5rj^^ |38| . We stress that, because we are dealing with model polymers, the solvent is 
not explicitly included in our simulation. The polymer-solvent interaction is effectively accounted for only through 
the Langevin thermostat and related friction term. 



4 



(a) 



(c) 




(d) 




FIG. 1. Typical configurations of ring polymer solutions, at pa^ = 0.025 (a) and pa'^ — 0.4 (b), and corresponding selected 
rings with trefoil-knotted (c and d, right) and unknotted topologies (c and d, left). 



Preparation of initial configurations - At each density, we consider monodispcrse solutions, consisting either of 
unknotted or trefoil-knotted rings. In the following we shall often use the standard shorthand notation UN and 3i to 
denote unknotted and trefoil- knotted topologies, respectively. 

The system initialization consists of placing a template unknotted or trefoil-knotted ring at the center of a cubic 
cell with periodic boundary conditions. The template configuration is replicated 4 times along each spatial direction. 
The linear dimension of the cubic cell is large enough to avoid overlap and linking between the template copies and 
yields an initial monomer density of about pa^ « 5 • 10~'^. Prior to the production runs, the cell is first evolved at 
constant pressure until the desired target density is reached and is subsequently equilibrated for a time span of about 
IO^tmd- At all considered densities, this time span exceeds the time required by one ring to diffuse over distances 
larger than its typical size. Production runs have a typical duration of 6 • 10^T7\/£). Typical equilibrated configurations 
at two different densities are shown in Fig. [1] 



As customary, the salient metric properties (shape and size) of a given ring are characterized through the gyration 
tensor, Q. The entries of this 3x3 matrix are given by: 



where a and /3 run over the three Cartesian components and rcM = SiLi spatial location of the ring center 

of mass. The non- negative eigenvalues of Q, ranked with decreasing magnitude, Ai, A2 and A3 correspond to the 
square length of the principal axes of the ring gyration ellipsoid. Accordingly, their relative magnitude conveniently 
captures the overall spatial anisotropy of the rings, while their sum yields the ring square radius of gyration: 



C. Ring shape and size 




(8) 
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(9) 



a=l 



where Tr is the trace operator. To measure the typical size and anisotropy of the ensemble of rings in solution, 
we compute the averages {Rg) and (Ai^i.2,3) where the brackets (...) denote averaging over all rings of all system 
snapshots at a given solution density. 
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(a) (b) (c) 



FIG. 2. Knotted portions in a non-trivially embedded trefoil knot (a), as identified by two different search approaclies. (b) 
Top-down approach: the red arc is the shortest continuously knotted portion. It can not be shortened without loosing the 
3i topology, (c) Bottom-up approach: the red is the shortest arc with 3i topology. Complementary arcs (in blue) in both 
approaches must be unknotted (see Ref. IssI], for details). 

D. Identifying the knotted portion of the rings 

Because of the chain connectivity constraint, the global topological state of the rings is preserved in the course of the 
MD evolution. One way to characterize the interplay of geometrical and topological entanglement is by establishing 
the magnitude of the typical length of the ring portion that accommodates the knot [H, Ha, SH . 

From a mathematical point of view, the problem of locating the knot in a closed chain is ambiguous [13, IH, \^ . 
The difficulties are essentially two [s^. First, to assess if the knot is contained in a given ring portion, which is an 
open arc, it is necessary to define a procedure for closing it by suitably bridging the arc ends. In fact, only after the 
arc is closed it is possible to characterize its topological state by using topological invariants. Viable closure procedure 
should introduce the least possible additional entanglement. This criterion is well satisfied by the minimally-interfering 
closure that some of us introduced recently [s^l and that will accordingly be used throughout our study. The second 
difficulty lies in the fact that the geometrical entanglement of a knotted ring can be so high that the identification of 
the region that accommodates the knot depends on the method used to search for it (s^ . Arguably, upper and lower 
bounds to the length of the region accommodating the knot are offered by the top-down and bottom-up searches 
illustrated in Fig. [51 Both methods were applied here and systematically yielded consistent results, see Supporting 
Fig. 1. 

Based on the robustness of the knot localization procedure for rings in solution, we shall limit considerations to the 
bottom-up search. This method returns the shortest arc having the same topology of the whole ring (a trefoil knot, 
in our case) while its complementary arc is unknotted, see Fig. 



E. Surface accessible area 

The rings surface accessible area was computed using the SAS routine of the GROMACS package [3, using a 
probe sphere of diameter equal to the bead diameter, a. The computed surface area includes contributions from inner 
cavities of the ring conformations that are large enough to accommodate the probe sphere. 



III. RESULTS AND DISCUSSION 
A. Geometric properties of the rings 

We first report on how various equilibrium metric properties of the rings depend on their topology and solution 
density. In a progression from global to finer aspects, we shall consider the rings size, shape, exposed surface area, 
the mean-square distance of points at increasing arclength separation and the degree of localization of the knot that 
they accommodate (in case of non-trivial topology). 

The rings consist of n = 216 beads of diameter a with associated Kuhn length, £k = 9(t. The considered monomer 
density of the solution, pa^, spans the interval [0.01, 0.4]. This range spans from a dilute situation to one where the 
inter-chain entanglement is significant and yet the solution is isotropic (46j . The lowest monomer density is much 
smaller than the one at which one expects significant ring overlap. The overlap density is estimated as Q, 

-~ (10) 



(i?2)3/2 



where {R'^) is the mean square radius of gyration of an isolated ring. Using the Zimm-Stockmayer estimate for the 
isolated ring size m, = ~ 162(7 yields p*a^ w 0.1, which is an order of magnitude larger than the minimal 

monomer density considered here. 
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Above this density, significant mutual entanglement of the chains is expected. For dense solutions of Zmear polymers, 
the intricacy of the melt @, @ and the resulting physical properties [ill are captured by the entanglement 
length, Lf,. Roughly speaking, this quantity corresponds to the typical chain arclength separation between two 
consecutive topological constraints (known as entanglements) arising from inter-chain wncrossability. It has been 
recently proposed |46| that the relationship tying Lg, in and the solution density, p is adequately captured by the 
following phenomenological expression: 

^ « (0.06 [patDr^/^ + (0.06 {pail))-\ (11) 

Using the previous expression, one has that the interval of p considered here corresponds to a wide range of Le for 
the equivalent system of linear chains. Indeed, for pa^ going from 10~^ to 0.4 we have that Le/^if spans from w 400 
to ?» 1. In the latter situation, the entanglement length is approximately equal to the Kuhn length, and hence the 
system is at the crossover from the loosely- to the tightly-entangled regimes (i.e. the chain is nearly straight between 
two consecutive entanglements [i^). Higher densities are not considered, because for pa^ > 0.5 the onset of the 
Onsager isotropic/nematic transition [i^ is expected to break the spatial isotropy of the solution. 





{Rl) [o^] 


(Ai> k'] 


(A2) [a^] 


(A3> k^] 


UN 
3i 


175.4 ±0.6 
120.8 ±0.7 


116.3 ±0.6 
77.9 ±0.7 


43.9 ±0.1 
30.6 ±0.1 


15.2 ±0.1 

12.3 ±0.1 


FJR 


168.7 ±0.4 


108.5 ±0.4 


42.1 ±0.2 


18.1 ±0.1 


GR 


161.6 ±0.7 


106.1 ±0.6 


39.7 ±0.1 


16.2 ±0.1 



TABLE I. Average value of the mean square radius of gyration and eigenvalues of the gyration tensor for unknotted and 
trefoil-knotted isolated rings. Analogous quantities are also shown for an equivalent freely-jointed ring (FJR) and Gaussian 
ring (GR) of 24 bonds. The FJR bond length and the GR root-mean square bond length are equal to Ik = 9a. 

Following the outline set at the beginning of this section we first discuss how overall ring metric properties depend 
on the monomer density, p. 

For reference, we provide in Table |T] the average values of the square radius of gyration, and the eigenvalues 

of the gyration tensor for isolated unknotted and trefoil-knotted rings. Analogous quantities are reported for an 
equivalent freely-jointed ring (FJR) and Gaussian ring (GR) of LcjiK = 24 bonds 0. The FJR data are obtained 
by a Markovian exploration of FJRs, while the GR data are obtained by stochastic molecular dynamics simulations. 
Table U conveys the effect that intra-chain constraints (excluded volume and fixed topology) have on the overall size 
and shape of rings at infinite dilution. Both types of constraints arc absent in the FJR and the GR. Notice, that the 
FJR and GR quantities arc typically within 10% of those of unknotted rings. The difference with analogous quantities 
for trefoil-knotted rings is substantially larger. These results reflect the fact that, for the considered values of and 
unknotted rings dominate the equilibrium ensemble of infinitely-thin (FJR) or Gaussian rings with unrestricted 
topology. Trefoil-knotted rings therefore possess a degree of entanglement that would be atypical in rings circularized 
in equilibrium. The presence of the non-trivial topological constraint causes such rings to be tighter and slightly 
more isotropic than unknotted ones. This result is consistent with the intuitive notion that a certain arclength of 
trefoil-knotted rings is "used up" in the non-trivial entanglement (in fact the minimal ropclcngth required to tie 
trefoil- knotted rings is larger than for unknotted ones). 





UN 


UN 


UN 


3i 


3i 


3i 


Monomer density, pa^ 


kRl) \A 


(Ai>/{A3> 


(A2>/(A3) 


kRl) k'l 


(Ai)/{A3) 


(A2>/(A3> 


0.010 


165.59 ±0.37 


7.45 ± 0.05 


2.77 ±0.01 


116.07 ±0.01 


6.118 ±0.001 


2.4151 ±0.0004 


0.025 


150.55 ± 0.39 


7.21 ±0.03 


2.63 ±0.02 


108.23 ± 0.08 


5.91 ±0.01 


2.335 ± 0.003 


0.050 


134.64 ± 0.06 


7.07 ± 0.02 


2.530 ± 0.003 


97.26 ±0.05 


5.80 ±0.01 


2.266 ± 0.006 


0.100 


115.63 ±0.09 


7.00 ± 0.03 


2.44 ± 0.07 


85.62 ±0.19 


5.95 ±0.02 


2.244 ± 0.004 


0.200 


98.64 ±0.10 


7.02 ± 0.02 


2.385 ± 0.006 


77.94 ± 0.22 


6.29 ±0.03 


2.260 ± 0.008 


0.400 


86.32 ±0.85 


7.46 ±0.13 


2.41 ±0.01 


70.11 ±0.06 


6.84 ±0.03 


2.318 ±0.008 



TABLE IL Mean square radius of gyration, (-Rg) and average shape of the simulated rings expressed as the ratios (Ai)/(A3) 
and (A2)/(A3), where (Aa=i,2,3) are the three average eigenvalues of the gyration tensor, Eq. [S] For very dilute conditions 
(pcr'^ — 0.010), the calculated quantities are close to the theoretical values for ideal semi-flexible rings (see Table 

Ring size. Table |lT] reports the same metric properties of Table HI but calculated for rings in solution. The data 
suggest that, for both topologies, the ring size decreases with the monomer density, p. Specifically, going from 
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pa^ = 10~^ to pa^ = 0.4 it is seen that i?^ decreases from the typical value of isolated rings down to about half of it. 
The decreasing trend of the average extension results from the non-concatenation constraint of the rings, consistent 
with previous numerical findings 0, 3 ■ It should also be noted that the monomer density attained by each ring is 
much smaller than the density of the entire solution. In fact, even in the densest case pa^ = 0.4 the individual ring 
density - defined as the number of ring monomers divided by the volume of the average gyration ellipsoid - is only 
« 0.05(7^'^ for unknots and « O.OGcr^"^ for trefoil-knotted rings. These findings suggest that an increase of solution 
density promotes extensive chain intermingling as opposed to tight compaction of the individual rings. 

Ring shape. A further indication that an increase of the solution density affects only weakly the metric properties of 
each ring is give n by looking at the typical ring shape. For this purpose, following previous studies of either closed or 
open chains |49l - [5l| . we calculate the average eigenvalues (Aq=i.2,3) of the gyration tensor. Eq. [8]which are shown in 
Table [III We note that the ratios (Ai)/(A3) and (A2)/(A3) fiuctuate by a few percent throughout the explored range 
of solution densities and remain close to the overall anisotropy expected for equivalent, isolated FJRs (see Table 

Exposed surface. To complete the statistical characterization of ring conformations we next focus on the average 
exposed surface area of each ring. This is an important geometric indicator that aptly complements the radius of 
gyration 0, @, 0j H!! in describing the overall degree of ring compactness attained in denser and denser solutions. 
We stress that, because the purpose is to characterize the properties of individual rings, the exposed surface per bead 
is calculated separately for each ring, i.e. without taking into account the burial effect due to the surrounding chains. 
Furthermore, the exposed surface is calculated by taking into account voids (cavities) possibly present in a ring (see 
Methods, Sec. UTE)) . 

The varying degree of exposure of the beads in a ring is illustrated in Fig. [SJj: the profile pertains to the unknotted 
ring conformation in panel (a) (picked at pa^ = 0.4), where the beads are color coded according to the degree of 
exposure. For reference, in panel (b) it is also shown the value of the exposed surface per bead in a perfectly straight 
chain configuration. For the case shown in panels (a-b), highly exposed beads nearly approach this reference value, 
while the most buried ones have an exposed surface equal to about 1/3 of the reference value. The distributions of 
the total exposed surface per ring are shown in panel (c): in particular, we notice that the ring exposed surface covers 
a moderately narrow interval. Nevertheless, the exposed area for trefoil-knotted rings is generally smaller than for 
the unknots (at the same monomer concentration, p), and the mean value of the exposed surface per bead (averaged 
over all beads of all rings in various configurations at fixed density) is a decreasing function of solution density. 

Finally, it is interesting to correlate the decrease of the surface area with the decrease of the ring radius of gyration 
observed for increasing density. This relationship is shown in panel (d) of Fig. [3] By taking into account the different 
scales and offsets of the two axes in the graph it is realized that unknotted and trefoil rings at the same monomer 
concentration differ appreciably by average size (see also Table |ll| though not by exposed surface. Secondly, the 
density-dependent decrease of the radius of gyration is not paralleled by an analogous decrease of the exposed surface. 
In fact, with respect to the diluted case, the latter diminishes by less than 5% up to pa^ = 0.2 and by less than 
10% up to pa^ = 0.4. We note that simple dimensional considerations would have suggested, instead, a proportional 
dependence of the surface accessible area on {B?g). The failure of the simple dimensional analysis suggests that, despite 
ring sizes that are, on average, smaller at higher densities, the corresponding configurations arc still sufficiently exposed 
to the solvent and maintain around each bead enough "free room" to accommodate the probe sphere used to measure 
the accessible surface. 

Geometry of ring portions. The results presented so far address properties of entire rings. We shall next examine 
how various ring portions, or arcs, are affected by intra-chain and inter-chain entanglement. We first report on how 
the mean square end-to-end distance, i?ee(^) of ^'^cs of contour length £, depends on the solution density and ring 
topology. Furthermore, for rings with non-trivial topology we shall identify the smallest arc accommodating the knot 
and examine how its contour length depends on the monomer density. 

The density-dependent behavior of i?ee(^) shown in Fig. 31 The data for knotted and unknotted rings are 
presented in two separate panels. For reference, in each panel we include the graph for the mean square end-to-end 
distance for a worm-like chain (WLC): 



2e 2£ 



2 

The analogous expression for a worm-like ring (WLR) is closely approximated by: 

/ \ -1 

1 



(12) 



V ^ee,WLCy'^) ^ee,WLC\^c ~ «j 



The latter expression is obtained by matching analytically the exact, small-£ (stiff) and large-^ (flexible) limiting 
behaviors for a ring polymer. 
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FIG. 3. (a) The shaded region outlines the accessible surface of the shown ring configuration. This is calculated using the 
SAS routine of the GROMACS package [3, SEl with a probe sphere of diameter equal to the bead diameter, a. Monomer 
beads are colored from red to cyan for increasing degree of their surface exposure, shown in panel (b) with the same coloring 
scheme, (c) Probability distribution functions of the total accessible surface per ring, for different solution densities and chain 
topologies, (d) Average accessible surface per ring as a function of the ring average size, for different solution densities and 
chain topologies. 




FIG. 4. Average square spatial distance i?^e between ring monomers as a function of their (normalized) contour length 
separation, £/Lc, along the chain: results for (a) unknotted and (b) trefoil-knotted ring polymers at increasing monomer 
densities. 
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Fig. |4] shows that for both unknotted and trefoil-knotted rings at ah densities, the WLC (and WLR) behavior is 
fohowed closely only for arclengths smaller than about £k (corresponding to i/L^ ~ 0.04). Beyond this arclength, 
the WLR behavior is followed well by unknotted rings in dilute solutions. Noticeable differences from the WLR 
trend occur for both unknotted and knotted rings at densities larger than p*a^ = 0.1. The deviations grow with; (i) 
increasing solution density at fixed topology and arclength separation; (ii) increasing arclength separation at fixed 
topology and density; (iii) changing ring topology from unknotted to trefoil-knotted at fixed density and arclength 
separation. In summary, departures from the WLR trend are seen upon increasing either the inter-chain (solution 
density) or the intra-chain (topological) entanglements. 

Properties of the knotted portion of a ring. We conclude the analysis of the equilibrium properties by examining 
how density influences the length of the ring portion accommodating the knot, ifenot- The analysis is motivated 
by recent studies which showed that the typical value of L^not is substantially larger for collapsed rings than for 
unconstrained or swollen ones (33.[36l[37|. In addition, for knotted rings subject to spherical confinement, it was seen 
that the increasing geometrical entanglement acquires a complex multiscale character so that very different values of 
Lknot are obtained upon locating the knotted portion using bottom-up or top-down approaches (see Methods, Sec. 
Ill pp . These properties are expected to significantly impact the physical behavior of the chains (see e.g. Ref. [sHi 
for a discussion of knot localization effects in a biological context, namely the ejection of tightly packed - and highly 
entangled - DNA from a viral capsid). 



Monomer density, pcr^ 


(Lknot) 


(Lknot)/ Lc 


0.010 


88.8 ±0.2 


0.41 ±0.0010 


0.025 


90.2 ±0.3 


0.42 ±0.0010 


0.050 


89.7 ±0.1 


0.42 ± 0.0005 


0.100 


85.9 ±0.1 


0.40 ± 0.0005 


0.200 


79.9 ±0.1 


0.37 ±0.0005 


0.400 


72.9 ±0.1 


0.34 ± 0.0005 



TABLE TIL Average size of the knotted portion of 3i ring polymers, (Lknot) as a function of monomer density, p. On average, 
the knot occupies a large fraction of the polymer chain. 

In the present context, the average value of L^not remains fairly constant at all densities, see Table Hill However, it 
should be noted that above the overlap density, p* , the probability distribution function for Lknot acquires a rounder 
and broader shape, see Fig. ??. The results are robust across different methods used for locating the knot (see 
Supporting Fig. 1 for a comparison of top-down and bottom- up approaches). In conclusion, the properties of knots 
in dense solutions differ qualitatively from the two most typical conditions (ring collapse or ring spatial confinement) 
leading to high monomer density. 
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1 

0.5 


0.2 0.4 0.6 0.8 1 
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FIG. 5. Probability distribution function of the normalized knot length, Lknot/ Lc, at different monomer concentrations. 



B. Dynamical properties 

The growing level of entanglement found in solutions of increasing density is expected to strongly impact the ring 
dynamics. By analogy with the case of linear polymers melts [ill [l3l l38l . l52j |. multiple dynamical regimes are 
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FIG. 6. Radius of gyration time correlation function CRgit) (Eq. I14|) for unknots (red symbols) and trefoils (green symbols), 
(a): Result for isolated rings (open symbols) and rings in very dilute conditions with pa'^ = 0.003 (filled symbols). The long 
time behavior of CRg{t) for 3i rings deviates from the Cr^ (t) for the unknots, and shows a longer exponential tail, in agreement 
with previous studies of isolated ring polymers on lattice [20l. [2^ [s^ . (b): Results for pa"^ — 0.01 (open symbols) and pa^ — 0.4 
(filled symbols). 




expected. For the latter system, it is known that at times longer than the slowest chain relaxation time, the motion 
of the chain center of mass follows the standard diffusion while at smaller times the kinetics is dominated by the slow 
chain reptation inside the tube created by inter-chain constraints [s^ . 
We accordingly monitor several observables: 

1. The typical time scale for ring diffusion in space; 

2. The time scales dominating the fluctuations in rings size and orientation; 

3. Finally, for rings with non-trivial topology, we examine how the knotted portion moves both in space and along 
the ring backbone. 

Some of the above mentioned quantities have been considered before in studies on melts of unknotted rings [^,| 
in lattice models of isolated knotted rings [20l . [23 . [s^ ] , and in polymer models of knot diffusion along DNA [5£ 
The present context therefore offers an opportunity to examine the impact of both solution density and ring topology 
on dynamics and compare the time-scales associated to the various phenomena. 

Time autocorrelation Junction for the radius of gyration. We start by considering the time autocorrelation function 
of an internal (i.e. independent of the ring absolute space position and orientation) quantity, namely the radius of 
gyration: 

{Rg{t)R,m-{R,)^ 

where the brackets, (...), denote the average over simulation time and over rings. The behavior of C/^^ (i) was considered 
in recent lattice studies of isolated., unconstrained knotted and unknotted rings [23 . [sij . These investigations showed 
that Cflg decays with a characteristic time, t^^, that is larger for knotted rings than for unknotted ones [U,!!^!- The 
effect reflects the enhanced self-hindrance of rings with non-trivial topology; consider, for instance, that the minimal 
number of crossings observed in any two-dimensional projection of an unknotted ring is zero, while it is three for 
trefoil-knotted rings (which hence must necessarily wind on themselves more than unknots). 

Investigating the decay properties of CRg{t) in the present context serves a twofold purpose. On the one hand, 
it can clarify whether the above-mentioned lattice results apply to ofF-latticc contexts, too. On the other hand, by 
analyzing the dependence of CR^{t) on solution density, it is possible to assess if, and to what extent, the increased 
intra-chain and inter-chain entanglement affect differently the size relaxation times of unknotted and knotted rings. 

According to this plan we first compute Cr^ from simulations of isolated knotted and unknotted rings. The results 
are shown in Fig. [S^ (open symbols) and fully support the lattice results regarding the slower relaxation of trefoils 
compared to unknots. For rings in solutions, at all considered densities {pa^ = 0.01 — 0.4), one is surprised to find 
that the Cr curves for trefoils and unknots display much smaller differences compared to the isolated case, sec Fig. 
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[6)3. Indeed, a behavior quantitatively similar to the isolated case is found only at densities lower than pcr^ = 0.003, 
as shown in Fig. [S^. The effect is quantitatively captured by calculating the effective relaxation time, tji^ , given by 
the integral of (t) , see Table IIVI The data in the table illustrate vividly that the density-dependent increase of 
the intra- and inter-chain entanglement obliterates differences in the average kinetic behavior of rings with different 
topological state. For both unknots and trefoils, Tfj^ has an order of magnitude increase over the considered range of 
solution density, p. Notice that upon increasing p, the value of tr^ for unknots overtakes the one of trefoils (see also 
the order of the red and green curves in panels a and b in Fig. (S]). 



Monomer density, pa'^ 


UN 


3i 

TRg 


UN 

'^diam 


3i 

'^diam 


UN 

TCM 


3i 

TCM 


Single ring 


(6.6 ±0.6) ■ 10^ 


(7.9 ±0.6) ■ 10^ 


(2.8 ±0.1) ■ 10^ 


(1.6 ±0.1) • 10^ 


3.1 ■ 10^ 


2.2 • 10^ 


0.003 


(7.0 ±0.6) ■ 10^ 


(7.0 ±0.6) ■ 10^ 


(2.7 ±0.1) • 10^ 


(1.6 ±0.1) • 10^ 


3.4 ■ 10^ 


2.4 • 10^ 


0.010 


(8.1 ±0.6) ■ 10^ 


(7.3 ±0.6) ■ 10^ 


(2.7 ±0.4) ■ 10^ 


(1.6 ±0.3) ■ 10^ 


4.8 ■ 10^ 


3.6 • 10^ 


0.025 


(1.1 ±0.1) ■ 10^ 


(7.8 ±0.6) ■ 10^ 


(2.8 ±0.4) • 10^ 


(1.6 ±0.3) • 10^ 


4.8 ■ 10^ 


3.6 • 10^ 


0.050 


(1.6 ±0.1) ■ 10^ 


(1.0 ±0.1) ■ 10^ 


(3.1 ±0.4) ■ 10^ 


(1.7 ±0.3) ■ 10^ 


6.0 ■ 10^ 


3.6 • 10^ 


0.100 


(2.7 ±0.1) ■ 10^ 


(1.5 ±0.1) ■ 10^ 


(3.9 ±0.4) • 10^ 


(2.2 ±0.3) • 10^ 


7.2 ■ 10^ 


4.8 • 10-' 


0.200 


(4.7 ±0.1) ■ 10^ 


(2.6 ±0.1) • 10^ 


(6.1 ±0.4) ■ 10=^ 


(3.3 ±0.3) ■ 10^ 


1.2 ■ 10* 


7.2 • 10^ 


0.400 


(1.25 ±0.01) • 10-* 


(7.3 ±0.1) ■ 10^ 


(1.59 ±0.04) • 10* 


(9.2 ±0.4) • 10^ 


3.9 ■ 10* 


2.2 • 10* 



TABLE IV. Correlation times r_Rg and Tdiam, calculated by numerical integration of the respective time correlation functions 
C'Rg{i) (Eq- llip and Cdiam{t) (Eq. I15|l . Numerical integration is limited to the time interval where correlation functions are 
> 10^^. TCM is defined as the typical time required by rings to transverse a region of linear size comparable to the ring average 
gyration radius, i.e. 5r^j\f (tca/) ~ (Rg) (see Eq. I16|) . All times are expressed in units of tmd- 

Reorientation time. We now turn to consider kinetic properties that do depend on the absolute orientation of rings 
in space. For linear polymers, it is customary to consider the autocorrelation function of the end-to-end vector [8| 
which, for closed chains, admits several generalizations Q. The one considered here is the time correlation of the ring 
diameter vector: 



{d{t) ■ djO)) 

Cdmm(i) = -^r- , (15) 

(lo'l ) 

where d is the vector joining two monomers that are "diametrically opposite" on the ring backbone, i.e. monomers with 
the largest possible chemical distance, Lc/2. Here the brackets (• • •) denote a multiple average: over the dynamical 
trajectory and over each pair of diametrically opposite beads for every ring in solution. Notice that Cdiam is sensitive 
to both changes in modulus of the diameter vector and to its absolute orientation. For this reason, Tdiam, the 
characteristic decay time of Cdiam{t) is customarily referred to as the reorientation time. 

Fig. [7] portrays Cdiam{t) for unknots (red symbols) and 3i knots (green symbols) at different monomer concentra- 
tions. Analogously to th^, we define the corresponding Tdiam as the numerical integral of Cdiam{t)- Final results are 
reported in Table HVl As for Cu^{t), the asymptotic decay time of Gdiam{t) depends on ring topology and is faster 
for trefoils than unknots. To the best of our knowledge this effect, that holds at all densities, was neither pointed out 
nor addressed before. Extending this analysis and considerations to rings with more complicated knot types could 
represent an interesting avenue for further investigations. 

Diffusion of the ring center of mass. Because trefoils are more compact than unknots they should diffuse faster 
in the solution, analogously to what they do in other types of media [stI . Issj . This intuitive expectation is indeed 
confirmed by inspecting Fig. [51 which portrays the mean square displacement of the ring center of mass, 5rQjy,j(t), for 
time lags of increasing duration, t: 

Sr^cMit) ^ {{rcAiit) - rcM{0)f)- (16) 

The displacement data given in Fig. [5] pertain to pa^ = 0.2; similar plots are obtained for different values of p. It 
is seen that over the time-span covered by the simulations, both unknotted and trefoil rings diffuse over distances 
that exceed by several orders of magnitude their average size. For traveled distances larger than the average ring 
size, the motion of the center of mass follows standard diffusion. The corresponding diffusion time tcm, defined as 
^''^"cKii'^CM) = (^g), is given in Table ITVl Notice, that the values of tcm are comparable to the values of Tdiam at the 
same solution density, consistently with the case of linear polymers [8| . 
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FIG. 7. Ring diameter time correlation function Cdiamit) (Eq- llSp for the unknots (red symbols) and trefoils (green symbols), 
(a): pa'^ — 0.01 (open symbols) and pa'^ = 0.025 (filled symbols), (b): pa^ — 0.1 (open symbols) and pa^ = 0.4 (filled symbols). 




FIG. 8. (a): Comparison between the time mean-square displacement of the rings center of mass {5r'^(^]^.j{t), Eq. I16|) for 
unknots (red symbols) and trefoils (green symbols) and the corresponding average square gyration radii (horizontal lines), for 
pa^ = 0.2. Within the considered time-span, rings diffuse over distances much larger than their typical sizes, (b): Diffusive 
motion of the center of mass of a randomly selected trefoil, for pa^ — 0.2. The red bar denotes the average ring size, {R^Y^"^ . 
Corresponding plots at different monomer densities look qualitatively similar. 



Backbone motion of the knot. Finally, we examine the dynamics in space and along the ring backbone of the knotted 
portion of the trefoils. This is a computationally-demanding task, because the numerically-costly identification of the 
knotted portion must be carried out for each ring of all sampled system snapshots. 

For simplicity, the instantaneous position of the knot on the ring is taken to coincide with the chemical coordinate 
of the midpoint of the knotted portion (defined with respect to the absolute monomer indexing of the ring set at 
the beginning of the simulation). The instantaneous knot midpoint position, Skit) is recorded at fixed time intervals, 
Ai = 12{]tmd. The incremental knot displacement between time snapshots n (= and n -|- 1, dk{n), is computed 
using the expression: 



dkin) 



Sk{t + At)~Sk{t) + -L, 



mod Lr> — 



2 



(17) 



where "mod" denotes the modulus operation. The previous expression ensures that the displacement is mapped in the 
[—Lc/2, Lc/2] interval. The value of At is suSiciently small that, at all considered densities, the occurrence of large 
"jumps" (say equal or larger than Lc/4 in modulus) is rare. This is illustrated in the inset of Fig. [5] (for pa^ = 0.1) 
which further highlights the Gaussian character of the distribution. 
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FIG. 9. (a) Time behavior of the knot midpoint position along the ring contour. Midpoint coordinates are defined with respect 
to a reference monomer chosen at the beginning of the simulation. On the right, four ring configurations sampled at different 
times along the trajectory, with the corresponding knotted portions marked in red. (b) Mean square displacement of the knot 
along the ring contour, Ss1{t) (Eq. I18|l . at density pcr^ = 0.1. The proportional dependence on time (red line passing through the 
origin) can be expressed as 2Dkt, where Dk is the one-dimensional diffusion coefficient: for the case shown, Dk = 0.21(7^ /ta/u 
(see Table [V] for a list of diffusion coefficients at all monomer concentrations). Inset: The Gaussian distribution (red curve) 
associated with the model one-dimensional diffusion process matches well the numerical distribution (black line) of incremental 
knot displacements, dk (Eq. I17p . 



The knot mean square displacement on the ring contour, Ssf.(t), is defined by 



(18) 



where the summation index runs over consecutive displacements, and the average (...) is taken (1) over the trajectory 
and (2) over the rings ensemble. Fig. [9] illustrates that a linear relationship holds between Ssf.{t) and time, indicating 
that the knot sliding on the ring backbone can be treated as a one-dimensional diffusion process. The corresponding 
diffusion coefficients, Dk, calculated at different solution densities are given in Table IVl which also reports the typical 
time, Tfc, required by the knot to diffuse by its average length on the ring backbone. Interestingly, Tk is almost 
unaffected by increasing monomer concentration. In addition, for a given solution density, exceeds all other 
characteristic times given in Table ITVl 



Monomer density, p 


Dk[a^ /tmd] 




0.010 


0.18 ±0.01 


2.1 ■ 10* 


0.025 


0.19 ±0.01 


2.0 • 10* 


0.050 


0.19 ±0.01 


2.0 ■ 10* 


0.100 


0.21 ±0.01 


1.8 ■ 10* 


0.200 


0.23 ±0.01 


1.6 ■ 10* 


0.400 


0.16 ±0.02 


2.3 ■ 10* 



TABLE V. Diffusion coefficients, Dk for the motion of the knotted portion along the ring contour length. The value of Dk 
corresponds to half the slope of the line that, passing through the origin, provides the best interpolation to Ssk{t) vs t. The 
typical diffusion (sliding) time, Tk, is defined as the time required to diffuse by a backbone distance equal to the typical knot 
contour length (« 40% of the ring contour length, see Table HIT)) . 

Motion of the knotted portion in space. The above observation indicates that the knot location on the ring backbone 
is practically quenched as a ring diffuses in space over a distance equal to its size. This suggests that the three- 
dimensional motion of the knotted portion should be akin to the one of a fixed, tagged, portion of the ring of 
comparable size. 

This observation was verified by considering the time dependence of the mean-square displacement of the knot 
center of mass, defined (analogously to Eq. as: 



SrlcAiit) = {{rkCMit) - rkCM{0)f). 



(19) 
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t [^mdI t ['^md] 

FIG. 10. (a) Comparison between the mean-square displacement of the center of mass of trefoils ((5r^M(i), Eq. [H ereen 
symbols) and the center of mass of their knotted portion (SrlcMit), Eq- 1191 brown symbols). The motion of the knotted 
portion is sub-diffusive at short times, and follows the chain global displacement at larger times, (b) Comparison between the 
mean-square displacement of the center of mass of unknots {SrQ]^j{t), Eq. 1161 green symbols) and the center of mass of a tagged 
portion (brown symbols) of linear size equal to the average knot length on 3i rings 0.4Lc, Table HIT)) . 

Fig. [TU](lcft panel) shows that the diffusive behavior of the observable VkCM sets in at times larger than the relaxation 
time of the whole chain (compare to Fig. [8|), while a subdiffusive behavior ~ with a « 0.7 is seen at smaller times. 
Notice that the same exponent, a w 0.7 was previously observed by looking at the stochastic dynamics of the center 
of mass of the knotted portion in isolated self-avoiding polygons on a cubic lattice [13]. Analogous plots at different 
monomer densities show the same subdiffusive behavior with an exponent a that decreases slightly with p (not shown). 
We have then calculated the mean square displacement of a randomly-picked ring portion spanning 40% of the contour 
length (the typical size of the knotted portion. Table of the unknotted ring. As shown in Fig. [TUb it was found 
that the tagged portion moves very compatibly with the motion of the knotted part, showing an analogous crossover 
from a diffusive to a subdiffusive behavior. This suggests that the underlying mechanism of the subdiffusive regime of 
the "quenched" knotted part is similar to the one that governs the subdiffusive behavior of any other (topologically 
trivial) tagged subregion of the ring. Note, however, that the effective length of the rings considered here is very small 
compared to the one of the lattice rings studied in [2J] , and no simple argument is presently available to predict if the 
a ~ 0.7 exponent is expected to be maintained for increasing contour lengths where the size of the knotted portion is 
supposed to be negligible with respect to the one of the whole ring. 

In conclusion, the above results provide a vivid picture of the key features of the dynamics of the knotted region. 
Specifically, they illustrate that the displacement of the knot center of mass follows "passively" the one of the whole 
ring because of the very long times required by the knot to slide along the ring. 



IV. SUMMARY AND CONCLUSIONS 

We reported on a systematic computational study of the equilibrium and dynamics of solutions of unconcatcnated 
ring polymers with different knot topology. Specifically, molecular dynamics simulations (with no explicit treatment 
of the solvent) at fixed-volume and constant temperature were carried out on bead-spring models of semiflexible 
unknotted and trefoil-knotted rings for several solution densities. This framework was used to explore the extent to 
which the interplay between topological constraints (knotting) and the geometrical self (intra-chain) and mutual (inter- 
chain) entanglements affect both the equilibrium and the dynamical properties of the rings. The study complements 
previous investigations of dense solutions of unlinked and unknotted ring polymers of various contour lengths. 

Regarding the equilibrium metric properties it is found that changes of the inter-chain or intra-chain entanglement 
operated by varying the solution density and/or ring topology, affect modestly the average ring size and shape 
compared to the infinitely-diluted case. Specifically, the root mean radius of gyration of both unknotted and trefoil 
rings at the highest solution densities (occupied volume fraction equal to 0.4) are at most about 40% smaller than 
in the unconstrained case. At all densities, trefoils are smaller (between 10% and 20% in linear size) and slightly 
more globular than unknots. Yet the average exposed surface of rings with different topology is practically the same 
and about constant at all densities. These results offer an interesting insight regarding the compactness of the rings 
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and, in particular, they indicate that the moderate decrease in ring size following the increasing intra- and inter-chain 
entanglement does not preclude the persistence of voids and cavities within the rings convex hull so that the exposed 
surface area is about the same as for unconstrained isolated rings. 

The weak dependence of the overall ring geometric features on solution density prompts the question of whether, 
for knotted rings, the metric properties of the knotted ring portion are weakly affected too. By knotted ring portion 
we refer to the shortest arc that accommodates the knot. This question is addressed by using the effective and 
transparent knot-location algorithm that was used recently to show that in isolated rings the length of the knotted 
region varies dramatically with increasing spherical confinement. At variance with the latter case, it is found that the 
length of the average size of the knotted ring portion is practically insensitive to density variations. 

While ring equilibrium metric properties are only weakly affected by variations of solution density and ring topology, 
the opposite holds for kinetic properties. The characteristic times of ring size relaxation, reorientation, and center of 
mass diffusion change by one order of magnitude across the considered density range. 

Further topology-dependent aspects of ring kinetics were highlighted by monitoring various kinetic observables of 
the ring knotted portion. Specifically we focused on how such region displaces in space and along the ring contour. 
The backbone motion is found to follow standard one-dimensional diffusion and the corresponding diffusion coefficient 
is smaller (much smaller at low concentrations) than the one of the ring center of mass. Consequently, for time scales 
over which a ring moves appreciably in space, its knotted portion remains practically quenched and, as wc verified, 
diffuses in space as any other equally-long tagged portion of the ring. 

In summary, the results presented here offer novel insights into the impact of intra-chain and inter-chain entangle- 
ment in solutions of ring polymers. The mild dependence of equilibrium metric properties on both the above-mentioned 
effects has no parallel with the behavior of dense systems of rings obtained by three-dimensional spatial (spherical) 
confinement. This suggests that an interesting novel avenue to address in future work would be to examine analogous 
effects for dense systems of rings obtained by two-and one-dimensional confinement (slabs and channels), which are 
increasingly adopted for advanced polymer micromanipulations. By converse, the sensitive dependence of various ring 
kinetic properties on solution density and ring topology suggests that both effects should be relevant to the studies 
of rheology and transport properties in semi-dilute solutions of circular polymers. 
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SUPPORTING FIG. 1. Three-dimensional histogram comparing top-down and bottom-up search schemes of the length, Lknot 
of the knotted portion of the chain (see Methods, Sec. IID): the histogram reports the number of times the pair of values 
given by (ife'^""'""'''', ^iTot ''°™") occur in the course of the simulation. It is evident that the two algorithms measure the same 
knotted portion in the majority of cases. The case shown corresponds to ring polymers at monomer concentration p = 0.4. 
Analogous results are found at the other concentrations. 



